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DIFFUSION ALGORITHMS AND DATA REDUCTION ROUTINE FOR ONSITE 
REAL-TIME LAUNCH PREDICTIONS FOR THE TRANSPORT OF 
DELTA-THOR EXHAUST EFFLUENTS 


SUMMARY 


Specialization of the National Aeronautics and Space Admin is tration/Marshall 
Space Flight Center (NASA/MSFC) Multilayer Diffusion Model for the Delta-Thor vehicle 
permits a programmable calculator (HP9820A) to be utilized at the launch site to make a 
forecast of the surface effects resulting from the exhaust effluents of this vehicle. This 
specialization of the general diffusion model limits the scope of analysis to the surface 
mixing layer and the ellipsoidal dispersion of the exhaust cloud (model 3) -- the resulting 
concentration and dosage prediction are limited to the earth’s surface. 

Because of the specialization, the only mandatory inputs to the Delta-Thor 
version of the NASA/MSFC Surface Layer Diffusion Model are the data points from a 
rawinsonde sounding and the altitude of the top of the surface mixing layer. The data 
reduction routine employs the calculator plotter to display the meteorological profile, the 
temporal history of the rise of the Delta-Thor exhaust cloud to the point of stabilization, 
and a ground- level delineation of the concentrations and dosage of hydrogen chloride 
(HCl) along the exhaust cloud path. In addition, the HCl concentration isopleths are 
plotted on a map of Kennedy Space Center and Cape Canaveral, Florida. 

The NASA/MSFC Surface Layer Diffusion Calculator Routine for the Delta-Thor is 
designed for both real-time launch-site predictions and as a guide in determining the 
inputs for the NASA/MSFC Multilayer Diffusion Computer Program. This calculator 
routine is similar to the calculator routine for the Titan III [ 1 ] , except for a more 
complex cloud rise routine and a change in source constants. 


I. INTRODUCTION 


The NASA/MSFC Multilayer Diffusion Algorithms have been specialized for the 
Delta-Thor vehicle into the NASA/MSFC Surface Layer Diffusion Calculator Program [2] 
to provide launch-site predictions of the rocket exhaust effluent transport. This calculator 
program has been designed for launches at Kennedy Space Center. 

Investigations to develop quantitative procedures for assessing the environmental 
effects from the aerospace release of effluents into the atmosphere were initiated in 1962 
at NASA’s Marshall Space Flight Center because it was recognized that universally 
accepted and adequately validated prediction techniques for the transport of exhaust 
effluents were not available. It rapidly became apparent that some uncertainty existed 
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concerning some of the fundamental aspects of this problem, such as the amount and 
composition of the rocket exhaust effluents, the deposition and scavenging of these 
effluents, and the dispersal and transport of these effluents in the atmosphere. The 
available atmospheric measurements to determine the reliability of the description of the 
rocket transport models in the atmosphere were sparse and of questionable accuracy in 
the 1960’s. In conjunction with the National Environmental Policy Act of 1969 and the 
April 23, 1972, guidelines of the Council on Environmental Quality, which required 
environmental assessments for aerospace operations, the need for implementing a NASA 
program for monitoring the transport of exhaust effluents from large rockets to 
determine the rehabiUty of the NASA/MSFC transport description was recognized. 

As a result of informal discussions between representatives of NASA 
Headquarters, MSEC, Langley Research Center (LaRC), and Kennedy Space Center 
(KSC), it became apparent that a NASA in-house rocket exhaust effluent prediction and 
measurement program was desirable, possible, and feasible. A joint solid rocket motor 
exhaust prediction (MSEC) and measurement (LaRC and KSC) program evolved in 1972 
utilizing the Titan and Delta launches as a source for empirical information that can be 
employed to more accurately predict the environmental effects from aerospace 
operations. 

The NASA/MSEC Multilayer Diffusion Model is an Eulerian transport model 
derived from the gradient diffusion theory [ 1 ] . The Cramer coefficients are employed to 
correlate the atmospheric parameters with the eddy diffusion coefficients in this model. 
Because of the comprehensive nature of this model, the large computation facility at 
MSEC was normally employed when the model was exercised in the initial research stages 
of this quantitative investigation. However, from the initial experience in launch support 
of the rocket exhaust effluent prediction and monitoring activities, it became apparent 
that real-time, onsite effluent transport predictions are desirable, if feasible. This 
feasibiUty implies that a mobile substitute for the large computational facility was needed 
to transform the description of the atmospheric kinematics and thermodynamics obtained 
from rawinsonde soundings and meteorological towers into a predicted effluent transport 
description. 

A mobile substitute for the large computer is a programmable calculator 
(HP9820A). Naturally, this means that the computer program, which requires a core 
storage of approximately 40K, must be reduced in size. This was achieved by taking the 
NASA/MSEC Multilayer Diffusion Model — which provides a transport description for 
any rocket exhaust effluents from any vehicle, anywhere in the troposphere, under any 
environmental conditions, and for different types of dispersive — and limiting the 
description’s appUcability to a specific effluent such as HCl, to a specific vehicle such as 
the Delta-Thor, to a specific layer in the troposphere such as the surface mixing layer, 
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and to spherical diffusion without deposition effects. The requirement for real-time 
predictions led to the incorporation of meteorological and buoyancy algorithms as an 
integral part of the NASA/MSFC Surface Layer Diffusion Program. In addition, user 
requirements dictated that the surface concentration footprint of the exhaust cloud be 
superimposed on a map of the Kennedy Space Center area. 

The modeling considerations utilized to achieve these results are discussed in 
Section II. The computational procedures and graphical results are discussed and provided 
in Section III. 


II. NASA/MSFC SURFACE LAYER DIFFUSION MODEL 


The spatial description, in terms of concentration and dosage, of the dispersive 
transport of effluents from a discrete source is afforded by the NASA/MSFC Multilayer 
Diffusion Model. Specifically, this application of the model is for the prediction of the 
distribution of the toxic effluents [hydrogen chloride (HCl), carbon monoxide (CO), and 
alundum (AI 2 O 3 )] associated with the rocket exhaust which are emitted during the 
launch of a Delta-Thor vehicle from Cape Canaveral, Florida. The prediction is used to 
assess the resulting environmental effects. The dispersive description afforded by this 
multilayer diffusion model is initiated at the point where the rocket exhaust 
cloud of effluents reaches thermodynamic equilibrium with the environment — cloud 
stabilization — and therefore depends strongly on the kinematic and thermodynamic 
profiles of the atmosphere, along with the chemical and thermodynamic composition of 
the exhaust cloud. 

Thus, the initial considerations in this section are given to the atmospheric 
description employed as the input to the cloud rise algorithms. The discussion of the 
stabilized cloud considers the algorithms for both the cloud rise and cloud dimensions 
utilized as the initial inputs to the NASA/MSFC Multilayer Diffusion Model. The final 
portion of this section provides the significant mathematical expressions that were 
employed to obtain the centerline concentrations and dosages and the concentration 
isopleths. 


A. Meteorological Profile 

The influence of the atmospheric conditions on the dispersive transport of rocket 
exhaust effluents is the basis for the format of the meteorological profile. To illustrate, 
consider the diffusion equation. The solution of the differential equation for the 
diffusion dictates that the diffusion problem be treated in two stages [2]. The nonlinear 
differential form of the diffusion equation is 


9 X(^ ,t) 
9t 


+ (v) 


V X (^.t) 


V- [K(T,t,p,T) • Vx(T,t)] 


( 1 ) 
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where 


X(r ,t) is the scalar concentration of the diffusing gas. 


<v> is the average wind velocity (temporal-spatial average). 

K(r,t,p,T) is the diagonal diffusion tensor which is a function of position, time, 
and the thermodynamic parameters, pressure and temperature. 


To solve this differential equation by the separation of variables, the necessary 
assumptions must be established to make this equation linear [3]. The normal technique 
— and the one used in the NASA/MSFC Multilayer Diffusion Model — to linearize the 
diffusion equation [equation (1)] is to restrict the model to a kinematic description by 
assuming that the diffusion coefficient (^) is a time averaged value and is 
thermodynamically independent. This implies the initial conditions occur when the rocket 
exhaust cloud achieves thermodynamic equilibrium with the atmosphere at cloud 
stabilization. Consequently, the effluent transport problem can be separated into two 
phases, the thermodynamic phase during cloud rise and the kinematic phase of diffusion. 


The atmospheric thermodynamic parameters (pressure, temperature, and density) 
govern the magnitude of the buoyant force on the exhaust cloud and thus dictate the 
height of cloud rise. Incorporating these atmospheric thermodynamic parameters into a 
suitable concise thermal description that will efficiently interface with the cloud rise 
algorithms requires the following considerations. Since the temperature (T) is a function 
of the pressure (p), as expressed in Poisson’s equation. 



( 2 ) 


where 


R is the universal gas constant. 

Cp is the specific heat at a constant pressure. 

R/Cp = 0.288. 

The concept of a potential temperature is introduced to reference the temperature to 
a specific pressure (1000 mb) and is defined as 


$ = 


T 



0.288 


(3) 
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The potential temperature can be shown to be a measure of the entropy 
(s= Cpln $ + constant); therefore, the vertical potential gradient (94>/9z) is a measure of 

the change in entropy [4] . Since in an adiabatic process there is not a change in entropy, 
the potential temperature gradient is zero (which corresponds to a straight vertical line on 
the meteorological profile). 

To achieve exhaust cloud stabilization with the atmosphere, an entropy balance 
must be achieved between the exhaust cloud and the atmosphere, which can be 
determined by utilizing the thermodynamic description afforded by the potential 
temperature profile. In the case of a hot rocket exhaust cloud, this balance results from 
both entrainment due to the turbulent mixing of this cloud and the exhaust cloud rising 
in the atmosphere to a region of higher entropy. If the potential temperature difference 
between the surface and a cloud height is negative or zero — what shall be defined here 
as an adiabatic condition — the entropy difference between the exhaust cloud and 
atmosphere will continue to increase and cloud stabilization will not occur. However, if 
the potential temperature gradient is positive — a stable condition — the entropy 
difference between the exhaust cloud and the atmosphere decreases as the exhaust cloud 
rises until equilibrium is obtained. Thus, the thermodynamic influences of the atmosphere 
on the hot rocket exhaust cloud during the initial transport stage where the exhaust 
cloud is rising to the point of equilibrium can be determined directly from the potential 
temperature profile. 

During the kinematic stage of the rocket exhaust effluent transport, the potential 
temperature profile can also be a guide to the top of the surface mixing layer in that the 
signature for the top of the surface mixing layer is characterized by a change in the 
meteorological conditions. Specifically, this signature can be a change in wind velocity, a 
change in thermal gradient, or both. The change in thermal gradient is normally a result 
of a temperature inversion or possibly an isothermal layer above the altitude of cloud 
stabilization — both of which result in a more positive potential temperature gradient. 

The NASA/MSFC Multilayer Diffusion Model depends strongly on the wind 
direction and speed profiles in the surface mixing layer (first 2 km) to determine the 
surface effects from the kinematics of the exhaust cloud. These profiles are required not 
only for the mean wind speed and direction at the cloud stabilization height, but also for 
the wind speed and directional shears over the surface mixing layer. 

The meteorological profile (wind speed, wind direction, temperature, and 
pressure) is normally obtained from a rawinsonde sounding of the atmosphere. To obtain 
the entropy profile required for these soundings, the dry bulb temperature and pressure 
are translated into the potential temperature in accordance with equation (3). 

If the levels of the sounding are referenced only in terms of pressure, rather than 
altitude (z) and pressure, the altitude of the level is computed using 
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(4) 


z 


Po 

29.3 (T + T„) In — 
« P 


where the subscript “o” refers to the next lower level that is referenced in the sounding 
to an altitude. 

Thus, the atmospheric thermodynamic parameters required for the cloud rise 
plume are the pressure and temperature which are utilized in obtaining the potential 
temperature profile. The atmospheric kinematic parameters required for the diffusion 
phase are the wind speed and direction. The delineation of these parameters as a function 
of altitude comprises the meteorological profile. 


B. Delta-Thor Exhaust Cloud-Rise Relations 

During a normal Delta-Thor launch, the burning of the rocket engine and motors 
results in the formation of a hot exhaust cloud of the rocket exhaust effluents. 
Subsequently, this exhaust cloud rises and entrains ambient air until thermodynamic 
equilibrium between the exhaust cloud and the ambient atmosphere is attained. The first 
stage of the Delta-Thor vehicle is comprised of a liquid-fueled engine and six solid rocket 
motors; therefore, the Delta-Thor has a moderate residence time on the launch pad — 
somewhere between that of the Titan [5] and the Saturn [6]. This characteristic results 
in a hybrid entrainment process that is a cross between spherical entrainment — an 
instantaneous source — and cylindrical entrainment — a continuous source. Since the 
propellant expenditure rate is approximately the same for the engine and the motors, the 
cloud rise relation for the Delta-Thor is, therefore, an average between the instantaneous 
and continuous models. In accordance with these assumptions, the cloud rise equations 
will be developed for the Delta-Thor vehicle that will afford similar results to those 
obtained by G. A. Briggs [7-]. The description developed here for the temporal 
delineation of the cloud rise during the thermodynamic stage is designed to afford the 
maximum height as well as the temporal delineation of the cloud rise. This temporal 
dehneation of the cloud rise is, in turn, coupled with the atmospheric kinematic profile to 
obtain the spatial history of the cloud rise in order to locate the point of the exhaust 
cloud stabihzation. 

The cloud rise formulation, as was pointed out previously, is dependent on the 
stability of the thermodynamic condition in the atmosphere. The potential temperature 
gradient (d<^/dz) is the index of the stabihty, which is given by 


R 

d$ _ / 1000 yP OT J_ ^ £ 
dz \ P / Cp T 
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where the potential temperature (^>) is defined in equation (3). The other variables in 
equation (5) are the dry bulb temperature (T,°K), the pressure (p, mb), and the altitude 
above the surface (z, m). The constants are the specific heat of air at constant pressure 
(Cp= 1.003 Joules/g°K), the gas constant (R = 0.289 Joules/g°K), and the gravitational 

acceleration (g, m/sec^). If 


A4> 

Az 


< 


0 


( 6 ) 


where A>J> and Az are the potential temperature difference and the altitude difference 
between the surface and the altitude of interest, then it is assumed that adiabatic rather 
than stable conditions exist. The cloud rise height (z) for a Delta-Thor for an adiabatic 
atmosphere is given by [8] 




-1% r- 


/rR 


3 Fj t^ /r 


R 


L2 7g"u \7g/J 


Va 


- + — 
T's Tg, 


- ’•rI- +- 


(7) 

whereas the cloud rise height (Zg) for a Delta-Thor for a stable atmosphere is given by 

[ 8 ] 
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The variables are the initial vertical cloud speed (Wq, m/sec), the buoyancy term 
(Fg, m'^/sec^ and Fg, m^/sec^), the stability term (s, sec'*), and the time after ignition 

(t,sec). The subscripts s and 1 denote solid or liquid. The constants are the entrainment 
coefficient (7s= 0.64 and 7g= 0.50) and the initial cloud radius (r^, m). 

The buoyancy term for the solid motor is 

3gQs 

F = ^ ^ , (9) 

^ ^'p Ps ^s 


and the buoyancy flux for the liquid engine is 


gQ. 


Fj = ^ 

^ c *r Pg Tg 


( 10 ) 


where 

Pg (g/m*) is the surface density. 

Tg (°K) is the surface temperature. 

Qg (cal) is the effective heat release. 

Og (cal/sec) is the rate of heat released. 


Based on a least squares curve fit of Titan IIIC data [2] , the heat release as a function of 
cloud height is 


Qs = 


0.39457 

1.32095 QgZR 


= 7.58557 X 


, . 0.39457 
10® z 


( 11 ) 


since 


Qg = 4.58812 X 10® cal/sec 


8 



Thus, the relationship between the buoyancy terms is 


F 


£ ~ 


1.009375 


-0.39457 

F,, z 


( 12 ) 


For the Delta-Thor, the initial cloud radius (rj^) is effectively zero, which permits 
equations (7) and (8) to be written as follows. 

1. Adbiatic exhaust cloud rise time: 


tA = i v/a (1 + 1 


259756 u z 


-0.605430 


(13) 


where 


» =y — = 


1.772533 X 10‘i‘ Pg Tg z 


3.605430 


(14) 


2. Stable case: 


ts = 


r 

sal 

-0.605430' 

arc cos 

+ arc cos 

2 J L 

- 0.629878 s a u z 


2v^ 


(15) 


where the stability term is 


g A4> _ 9.8 AcE) 

^ " Tg aT " t7 aT 


(16) 


The potential temperature gradient is taken directly from the potential temperature 
profile at launch time. 
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The exhaust cloud rise algorithm for the Delta-Thor [equations (13) and (15)] 
affords only a temporal description of the ascent of the exhaust cloud to the point of 
cloud stabilization. Examination of the adiabatic cloud rise relation [equation (13)] 
reveals that a time for cloud stabilization is not defined, and this may be disturbing to 
some. However, two points should be remembered: 

1. When there is an adiabatic lapse rate in the atmosphere, this lapse rate exists 
in a relatively thin layer near the ground, with a stable layer above it. 

2. Since, as was pointed out previously, an entropy balance between the exhaust 
cloud and its local environment is required, cloud stabilization only can occur where 
there is a positive potential temperature difference between an altitude and the surface; 
that is, in a stable layer. The altitude range in kilometers of the stable cloud rise 
algorithm [equation (15)] is between the surface and cloud stabilization; i.e., 

-0.277 r -O.O 5 I 

0 < z < (pgTgS^) j^l + 4u ^ (pg Tg) J km (17) 

Beyond these limits, the periodic nature of the algorithm results in variations outside 
these hmits. Since the stability term (s) is a function of altitude [equation (16)], the 
maximum cloud rise height requires some form of an iterative solution to obtain this 
height. 


To summarize, the temporal cloud rise delineation that is obtained from the cloud 
rise algorithms is interfaced with the meteorological kinematics to predict the cloud path 
before stabilization. It should be noted that the cloud rise algorithms given in this part 
are vahd only for the Delta-Thor since the source strength utilized is for the Delta-Thor. 
In addition, these cloud rise relations are a hybrid cross between the instantaneous and 
continuous cloud rise relations. 


C. Delta-Thor Version of NASA/MSFC Multilayer Diffusion Model 

By specialization of the NASA/MSFC Multilayer Diffusion Model for the 
Delta-Thor exhaust effluents, the prediction for the ground-level concentration isopleths 
can be obtained from a small programmable desk calculator (HP9820) in real time. The 
modeling approach employed is as follows. 

The general differential equation for kinematic diffusion [equation (1)] can be 
linearized by assuming that the meteorological profile represents the average atmospheric 
conditions over the area of interest and solved by separation of variables for the spatial 
distribution of the concentration and dosage resulting from the launch of the Delta-Thor. 
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The temporal variations are accounted for by the standard deviation terms for the 
elevation and azimuth. Because of the complexity of the resulting formulation, a 
generalized model will first be presented as an introduction to the actual algorithm [2] . 

The generalized concentration model for a nearly instantaneous source is 
expressed as the product of five modular terms: 


Concentration = {Peak Concentration Term} X {Along-Wind Term} 

X { Lateral Term } X {Vertical Term} X {Depletion Term} , 


whereas the generalized dosage model for a nearly instantaneous source is defined by the 
product of four modular terms: 


Dosage = {Peak Dosage Terms} X {Lateral Term} X (Vertical Term} 
X { Depletion Term} 


Thus, the mathematical description for the concentration and dosage models permits 
flexibility in application to various sources and for changing atmospheric parameters 
while always maintaining a rigorous mass balance. 

Two obvious differences exist. First, the peak concentration term refers to the 
concentration at the point x, y = 0, z = H (where x is the wind direction and H is any 
height) and is defined by the expression 


Peak Concentration 


Q, 


m 


3/ 

(27 t ) Oy 


(18) 


where Qj^ is the mass source strength and is the standard deviation of the 
concentration distribution in the ith direction. The peak dosage term is given by 

Peak Dosage = — — , (19) 

2 7T U Gy 
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where u is the mean wind speed over the layer. The second difference between these 
models is that the concentration contains a modular along-wind term to account for 
downstream temporal effects not considered in the dosage model. The along-wind term 
affords an exponential decay in concentration as a function of cloud transit time, 
concentration distribution, and the mean wind speed. 

The lateral term (which is common to both models) is another exponential decay 
term and is a function of the Gaussian spreading rate and the distance laterally from the 
mean wind azimuth. The vertical term (again common to both models) is a rather 
complex decay function since it contains a multiple reflection term for the point source 
which stops the vertical cloud development at the top of the mixing layer and eventually 
changes the form of the vertical concentration distribution from Gaussian to rectangular. 
The last modular in both models is the depletion term. This term accounts for the loss of 
material by simple decay processes, precipitation scavenging, or gravitational settling. The 
depletion term will be neglected in the Delta-Thor version of the model. 

The meteorological profile is utihzed in layering the atmosphere in accordance 
with homogeneous kinematic and thermodynamic properties — hence the name 
“multilayer diffusion model.” The specialization of the general NASA/MSFC Multilayer 
Diffusion Model has limited the number of layers for consideration to just the surface 
mixing layer. In addition, it is assumed that the source has a spherical shape with 
ellipsoidal expansion (model 3 of the NASA/MSFC Multilayer Diffusion Model [ 1 ] ). 

Thus, the dosage algorithm is 


D {x, y, zg < z < 


<m 


2tt Oy u 


exp I 




exp 


-(H-z)^ 


+ exp 
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2^z' 
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( 20 ) 
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where corresponds to the source strength or total mass of material in the layer, H is 

the height of the centroid of the stabilized cloud, and subscripts T and B stand for the 
top and bottom. 

By restricting the dosage mapping to the surface and defining the bottom of the 
layer as the surface, equation (20) simplifies to 


D {x, y, o } = 


In Gy a^u 


exp 


(#)!•■'© 


n 

-s 

i=l 


exp 


- (2i zj + H)2 


+ exp 


- (2i zj - H) 


( 21 ) 


where is the altitude of the top of tlie surface mixing layer and n is such that the first 

exponential in the summation is greater than 100. This is the specialized dosage algorithm 
that is used. 

The source strength (Q, ppm) for the HCl from the Delta-Thor is 


Q = Wg • 


10^ f 


HCl 


/ 22.4 \/ \ /1013.2\ 

[mhCiJv73.16J[ p, J 


( 22 ) 


where 

n '?94S7 

Wg is the fuel expenditure rate of the solid motors (7.983386 X 10® z g). 

fjjQ is the factor of HCl in propellant (0.2083). 

Mhci is the molecular weight of HCl (36.5). 

The source strength is 

Q = 3.78538 X 10* — X , (23) 

Ps 
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where the surface temperature (Tg, °K) and the surface pressure (Pg, mb) are used, since 
only the concentrations and dosages at the surface are of interest. 

The standard deviation of the crosswind distribution (ffy), the standard deviation 
of the vertical distribution (a^), and the mean wind speed (u) are defined as follows. 

1 . The standard deviation of the crosswind dosage distribution is defined by 



(24) 


where o^(t) corresponds to the mean layer standard deviation of the wind azimuth for 

the cloud stabilization time (t). The difference in the wind direction (A0', radians) is 
taken between the surface and the top of the surface mixing layer in accordance with 


Ad' — (d'j’ — 



(25) 


where 0-p and are the mean wind direction in degrees at the top and at the base of 

the layer, respectively. This is the wind shear. If the diffusion coefficient is again assumed 
to be one (a = 1 ), equation (24) becomes 



From this relation, one can determine the importance of the wind shear in determining 
the crosswind distribution of the effluent. In the surface layer 



(m+l)(zj-z) (z) 


(27) 
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where the standard deviation of the wind azimuth angle height Zj^ and for 

the cloud stabilization time t is 


‘^XR{r}= aAR{ro> 




( 28 ) 


Here <^ar{’’o} standard deviation of the wind azimuth angle at height zr and for 

the reference time period (Tq), and the power-law exponent (m) for the vertical profile of 
the standard deviation of the wind azimuth angle in the surface layer is 



Then, 


^'atM ~ <^AT (’’o) 



( 30 ) 


where cfat('^o} standard deviation of the wind azimuth angle at the top of the 

surface layer for the reference time period. 

The standard deviation of the wind azimuth angle cannot be determined at the 
top of the surface mixing layer because a temporal history of the wind azimuth cannot 
be obtained. This is normally calculated using the power-law relation 


{’’o} 



( 31 ) 
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where 


- ■ “-(S)/-© • 

where u is the wind speed at the top of the surface mixing layer. 
The crosswind virtual distance is 


*^yo 


(33) 


where Oy^ is the standard deviation of the lateral source dimension, which is 


Uyo = 7 H = 0.64 H 


(34) 


where H(m) is the height of the exhaust cloud stabilization. 

2. The standard deviation of the vertical dosage distribution is defined by the 
expression 


" ‘^E^Rz 




(35) 


where 

0 £ describes the mean standard deviation of the wind elevation angle, 
gives the vertical virtual distance, 
accounts for vertical diffusion. 

X £2 is the distance over which rectilinear vertical expansion occurs downwind 
from an ideal point source. 
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In this specialization, it is assumed that the vertical diffusion coefficient is one 
(/5=1), which permits rewriting equation (35) as 


Oz = o'e (x + Xz> 


(36) 


where 


u'e = o^A • 

Thus, is obtained from equation (32), since (in general) the value for ot^ cannot be 
obtained without special instrumentation. 

The vertical distance Xz is given by the expression 


X 


z 



(38) 


where Uzq = Oy^ is the standard deviation of the vertical dosage distribution at the 
distance from the source where the measurement is made in the surface mixing layer. 

3. The mean speed of cloud transport (u) in the surface layer is defined in 
accordance with the power law 



where u^ is the, mean wind speed measured at the reference height z^ and the 
power-law exponent (p) for the wind speed profile in the surface layer is described by 


P = log 



(40) 
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Here, uj corresponds to the mean wind speed at the top of the surface layer (zj). Thus, 
in the surface layer, the mean cloud transport speed ( u ) is 


Un 

J tP dz , (41) 

(Zt-Zr)ZrP 2r 


which reduces to 


(Zj-Zr) (zr)P (1+p) 


(42) 


The concentration (x) follows directly from the results for the dosage (D) 
algorithm given in equation (21). The average concentration then is just 


x{x,y,z} = D{x,y,z} 


(43) 


where the standard deviation of the along-wind concentration distribution in the 
layer is 




V2 


(44) 


and the along-wind cloud length (L{x}) for a point source in the layer at the distance x 
from the source is 


L{x} = 


0.28 (A u) (x) 
u 


Here, Au is the vertical wind speed shear in the layer and is defined as 


(45) 


A U = Uj - Ur 


(46) 
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and is the standard deviation of the along-wind source dimension in the layer at the 

point of cloud stabilization. The above equation for L{x} is based on the theoretical and 
empirical results reported by Tyldesley and Wallington [9], who analyzed ground-level 
concentration measurements made at distances of 5 to 120 km downwind from 
instantaneous line-source releases. 

In summary,' it should be pointed out that the standard deviations of the vertical, 
crosswind, and along-wind terms represent the cloud dimensions (Lj); that is. 


Li = 


4.3 a 


1 


(47) 


The factor 4.3 represents the 97-percent confidence level of a normal distribution. Hence, 
the initial source dimension is translated into the standard deviation initially for 
modeling. The standard deviations a^, Oy, give the cloud size during the diffusion 
process. 


III. NASA/MSFC SURFACE LAYER DIFFUSION CALCULATOR PROGRAM 


The NASA/MSFC Multilayer Diffusion Calculator Program has been specialized to 
provide the NASA/MSFC Surface Layer Diffusion Computer Program in accordance with 
the discussion given in Section II. The specialization amounts to restricting the 
applicability to the surface prediction of the transport of HCl for the Delta-Thor. In 
addition to these modifications, the NASA/MSFC Surface Layer Diffusion Program also • 
incorporates algorithms to calculate a meteorological profile and cloud rise. 

Data requirements are a rawinsonde sounding from the surface to approximately 2 
km (or about 7000 ft) and the standard deviation of the wind azimuth at the surface. 
This program utilizes this information to generate graphs for; (1) meteorological profile, 
(2) the temporal history of the exhaust cloud rise to the point of stabilization, (3) the 
surface centerline concentration and dosage of HCl along the cloud path, and (4) the 
surface HCl isopleth along with this path. The operations of the NASA/MSFC Surface 
Layer Diffusion Program are straightforward since inputs are called through calculator 
display. The overall block diagram, showing how the routines are interfaced, is given in 
Figure 1. 

The initial discussion in this section affords a description of the operation of the 
program along with its inputs and outputs. The second part gives detailed instructions for 
operating the NASA/MSFC Surface Layer Diffusion Program. These instructions are 
summarized in Appendix A. The last part presents the operational sophistication for 
processing time reduction. 
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Figure 1 . Block diagram of NASA/MSFC Surface Layer Diffusion Program. 


A. Overviews of the Computations, Inputs, and Outputs 
Performed in Each Routine 

The NASA/MSFC Surface Layer Diffusion Program is divided into three parts to 
meet storage limitations of the HP9820 calculator. The first part is the MET routine, the 
second part is the cloud rise routine, and the third part is the diffusion routine (Fig. 1). 
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During operations, the MET routine, which generates the inputs for the diffusion model, 
is erased after cloud rise computations with the loading of the diffusion routine. The 
meteorological and cloud rise data in storage are not affected during this transition. 

Both of these routines can be subdivided into sections where a computational or 
plotting procedure is carried out. The MET routine performs four processes — two are 
involved with raw data processing, and the other two with plot generation. These sections 
are the sort, convert, meteorological profile, and cloud rise. Similarly, the diffusion 
routine contains four sections of which three sections are involved in generating the two 
diffusion mappings and the other section contains the basic algorithms for the 
concentration and dosage calculations. These are the centerline dosage and concentration 
section, cloud path section, isopleth section, and the FA subroutine. The FA subroutine 
differs from the other sections in that this subroutine is utilized more than once in part 
of the program. 

The first subroutine of the MET routine is the sort subroutine, which takes the 
rawinsonde sounding (Table 1), sorts out the data sets (pressure, temperature, altitude, 
and wind velocity) by pressure, and places these sets in descending order into storage. 
Termination factors for ending processing loops are derived at the completion of data 
entry for use in the MET routine. 


TABLE 1. INPUT DATA FOR THE SORT ROUTINE 


TEST NBR 0475 
RAWINSONDE RUN 
CAPE CANAVERAL MTA. FLA. 
12/20/73 1545Z 
ASCENT NBR 0092 


ALT FT. 

WDIR 

WKTS 

TEMP 

DEW PT 

PRESS 

RM 

AS HUM 

DEN 

IR 

VS 

WS 

16 

200 

12 

17.0 

10.8 

1022.0 

66 

9.63 

1220.8 

330 

668 


1000 

213 

19 

13.5 

5.3 

986.5 

57 

6.74 

1194.8 

307 

659 

12 

2000 

236 

15 

11.6 

2.5 

951.4 

58 

5.96 

1160.2 

295 

657 

8 

3000 

211 

12 

15.1 

-16-6 

917.6 

9 

1.25 

1108.0 

254 

661 

6 

4000 

270 

10 

14.5 

-18.1 

885.1 

8 

1.10 

1071.2 

245 

660 

2 

5000 

259 

5 

13.7 

-19.4 

853.7 

8 

0.99 

1035.8 

236 

660 

4 

6000 

232 

3 

12.1 

-20.4 

823.2 

8 

0.91 

1004.6 

229 

658 

2 

7000 

237 

4 

9.9 

-22.4 

793.6 

8 

0.77 

975.9 

222 

655 

1 


MANDATORY LEVELS 


ALT FT. 

WDIR 

WKTS 

TEMP 

DEWPT 

PRESS 

RM 

626 

223 

15 

14.8 

7.4 

1000.0 

61 

2043 

237 

15 

11.8 

1.5 

950.0 

56 

3539 

213 

11 

14.7 

-17.4 

900.0 

9 

3120 

256 

5 

13.6 

-19.5 

850.0 

8 

6751 

233 

4 

10.5 

-21.8 

800.0 

8 


1 SIGNIFICANT LEVELS ] 

TEMP 

DEWPT 

PRESS 

17.0 

10.8 

1022.0 

13.4 

5.2 

986.0 

11.2 

4.7 

955.0 

9.0 

-23.5 

782.0 
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ALTITUDE (m) 


The second subroutine is the convert subroutine, which takes the data from the 
sort subroutine and converts the units of the parameters into units that are compatible 
with the calculations. Unit conversions are feet to meters (altitude), knots to 
meter/seconds (wind speed), and dry bulb temperature to potential temperature. If the 
altitude for a data set was not given in the rawinsonde sounding, the altitude is calculated 
using equation (4). The dry bulb temperature is plotted in this section (Fig. 2), and 
options exist at the end of this subroutine to record the processed meteorological data. 
The data from a formerly processed sounding can be reanalyzed by introducing the data 
at the end of this subroutine. 



WIND DIRECTION (deg) 


180 220 260 300 


Figure 2. Meteorological profile. 

The third subroutine, the meteorological profile subroutine, plots out wind 
direction, wind speed, and potential temperature (Fig. 2). This routine is strictly 
peripheral to the diffusion calculations since it does not contain any data processing. 

The next routine is the cloud rise routine, which affords the temporary history of 
the rise of the exhaust cloud during the initial 15 minutes of transport (Fig. 3). The 
kinematic and geometric parameters associated with the stabilization of the exhaust cloud 
are stored for later operations. To prevent the necessity of repeating this data processing 
in the first two routines, the stored results can be recorded (RFC “DA”, R102). 
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TEMPORAL ASCENT OF A DELTA-THOR EXHAUST CLOUD 

RAWINSONDE SOUNDING 12/20/73 1545Z CAPE CANAVERAL, FLA. 
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Figure 3. Temporal ascent history of exhaust effluents. 

The first subroutine of the diffusion routine is the concentration and dosage 
subroutine. This subroutine generates concentration and dosage levels over a 20-km 
distance, initiating at the point of cloud stabilization (Fig. 4). Maximum concentration 
and the height of the surface mixing layer are printed out. It is often helpful to record 
the results in the calculator at this stage. The second subroutine is the cloud path 
subroutine which predicts and plots the exhaust cloud’s transit path on a graphical map 
of the launch area (Fig. 5). This same graph also incorporates the predicted HCl 
concentration isopleths. Indices with time relations are generated to show the location of 
exhaust cloud stabilization, the point of maximum concentration, and the cloud’s arrival 
time at some arbitrary down range distance. The third subroutine, the isopleth 
subroutine, is used strictly for plotting, since it illustrates where contours of constant 
concentration of HCl reside in the launch area. This subroutine recycles without 
termination, which permits several isopleths to be plotted. Termination is not necessary 
since this is the last operational portion of the program. 

The FA subroutine in the diffusion routine contains portions of the specialized 
version of model 3 from the NASA/MSFC Multilayer Diffusion Models. Factors produced 
by this subroutine are combined by other subroutines to form concentrations and dosage 
values. 
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Figure 4. Centerline concentrations and dosages for hydrogen chloride. 

Data requirements for the surface layer diffusion calculator program fall into two 
categories: (1) the parameters governing the exhaust cloud source as defined by the 
vehicle and (2) the meteorological parameters of the surface layer obtained from a 
rawinsonde sounding and surface measurements. 

The exhaust cloud source data for the Delta-Thor are programmed into the 
command lines of the program as a portion of a constant to its respective equation. In 
Delta-Thor case studies, reentry of these parameters is not necessary. Altering source data 
requires the reprogramming of command lines in accordance with the discussion in 
Section II. 

The four parameters that make up the source data are: (1) the heat released by 
the Delta-Thor exhaust [Qjj, equation (11)] and utilized in the buoyancy terms [F, 

equations (9) and (10)] (this is located in line 3 of the cloud rise routine); (2) the 
fractional amount of material expelled into the surface layer (f) and required for 
calculating the mass source strength [Qj|j, equation (23)] (this is located in line 0 of the 

diffusion routine); (3) the entrainment coefficient [7, equations (7) and (8)] utilized in 
the cloud rise relations and in the calculations for the standards in wind shear Oy^, 

a^Q, equation (47)] (this is located in line 3 of the cloud rise routine and in line 1 of the 
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diffusion routine); and (4) the molecular weight of the toxin of interest (M) and required 
for calculating Qj^ (this is located in line 1 of the cloud routine). Respective values for 

Qjj, f, 7, and M for the Delta-Thor in HCl measurements are 7.585571 x 10® ,0.2083, 

0.57, and 36.5 [8]. 

Meteorological data are entered into the HP9820 computer through the keyboard 
on display calls from the program. The program offers data card loading and recording to 
facihtate operations in case study reproductions and in forming data references. The five 
meteorological profile parameters utilized by the surface layer diffusion routine are 
altitude with its respective wind direction, wind speed, temperature, and pressure. Units 
for these parameters can be either English or metric, since the necessary conversion 
routines are contained within the program. The only data which impact the diffusion 
prediction are within the surface layer, which is normally contained in the first 2.1 km of 
the atmosphere. These five meteorological parameters constitute one data set. The 
maximum number of data sets that storage can handle is 19, although at least 15 are 
required to prevent processing overflows. Generally, from 15 to 18 sets of data can be 
derived from one rawinsonde recording. In cases where less exist, the program is designed 
to be recycled without a data entry until 15 levels of information are in the machine. 
The date of the rawinsonde data is requested by the program to be used for referencing 
the outputs. The input occurs following meteorology profile parameter plots. There are 
no data inputs required for the cloud rise routine. 

The inputs for the diffusion routine correspond to wind direction, wind speed, 
altitude, and standard deviation of wind azimuth at ground level. These parameters are 
selected from the meteorology profile to indicate where the top of the surface layer 
mixing is. In general, this should be at a level that was entered into the calculator. The 
last input appears in the center of the isopleth routine. This requires a concentration 
value with which its corresponding isopleth is to be plotted. 

The output generated by the NASA/MSFC Surface Layer Diffusion Calculator 
Program is designed to afford both a spatial and temporal prediction for the footprint 
from the transport of the solid rocket effluents. In addition, values that are important in 
the analysis of the resulting predictions are printed on these graphs. To reduce the data 
reduction time, separate programs exist to format each graph with respect to setting up 
the proper coordinate system and the proper labelings (in normal operation, it has been 
found that reproducing these from the master is a time saving operation). 

The meteorological profile represents the atmospheric kinematics and 
thermodynamics that are relevant to the dispersion of the rocket exhaust effluents. The 
potential temperature is given in this profile because it is the primary determining 
parameter in predicting the altitude at which the rocket exhaust cloud will stabilize. The 
selection of whether the adiabatic or the stabilized cloud rise relation is utilized in the 
calculation of the cloud stabilization height depends on the adiabatic lapse rate relative to 
the surface. The adiabatic potential temperature regime, the stable potential temperature 
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regime, and the stable potential temperature are defined on the meteorological profile. 
The dynamic range for wind speed and temperature (potential and dry bulb) profiles is 
from 0 to 36.5 m/sec and from -4®C to 36.5^C, respectively. Wind direction scaling varies 
in increments of 20 deg to actual conditions with a maximum range of 0 to 360 deg. The 
program calls for a stop after each profile is plotted so that the color of the pen can be 
changed, thereby making each profile easily distinguishable from the other profiles. In 
addition to the profiles, the surface conditions (such as temperature, pressure, and 
density), the date, and the time of the rawinsonde soundings are printed on the 
meteorological profile. 

The cloud rise graph affords the temporal history over the initial 15 minutes after 
rocket ignition. The type of atmosphere (adiabatic or stable) is indicated in the 
delineation of the cloud rise. The exact values for the altitude and transport time to 
cloud stabilization are printed out to aid the investigator. 

The third graph routine outputs the dosage and concentrations of HCl along the 
centerline of the exhaust cloud path. This delineation ranges over the first 20 km after 
cloud stabilization. The height of the surface mixing layer and the maximum 
concentration of the HCl are printed out. The concentration and dosages of other rocket 
exhaust constituents can be obtained by multiplying by the constants of proportionality, 
which are: 

2.93 for CO (ppm) 

2.89 for AI 2 O 3 (mg/m^) 

The maximum values of the concentration and dosage that can be drawn are 9.5 ppm 
and 950 ppm-sec. 

The final output is a graph of the HCl isopleths of constant concentration. This 
routine plots the actual cloud path showing where the cloud stabilization and the 
maximum concentration occur. In addition, an arbitrary point is marked with the time of 
cloud passage to enable the user to determine the passage time at any point (note that it 
is assumed that the cloud has a constant velocity after cloud stabilization). Numerical 
printouts give the time of cloud passage, the distance from the launch site, and the 
Cartesian angle to the point of maximum concentration. 

The details of the operation of this calculator program in the next subsection will 
help illuminate points that are not clear. 
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B. Operating Instructions for the NASA/MSFC Surface Layer 
Diffusion Calculator Program 


A detaUed set of operating instructions for the NASA/MSFC Surface Layer 
Diffusion Program is provided in this subsection. While this set of instructions is helpful 
when initially running the program, the summary of operating instructions given in 
Appendix A should be sufficient in normal operations. The NASA/MSFC Surface Layer 
Diffusion Program is loaded into the machine (HP9820) in three parts, with the second 
and third parts depending on the preceding part. Each part has two sections, each 
affording one of the four graphs. The normal operation of this program can be 
summarized as loading the MET routine and then running the program. Data will be 
called for in the display. The user must terminate the data entry with a GO TO 12. 
Otherwise, all other instructions are called for until the meteorological profile graph is 
completed. Then the cloud rise routine is loaded. Loading of the third part (the diffusion 
routine) requires a bit more care, since the FA subroutine must also be entered. Once 
this is completed, the remainder of the program is straightforward. 

In the following instructions, the names of calculator keys are underlined and the 
displays are both underlined and set off by quotation marks. 

PARTI. MET ROUTINE 


SECTION I. SORT ROUTINE 

A. Loading MET Routine 

1. The MET routine consists of two cards, three sides, numerically listed for 

loading. 

2. Press ERASE MEMORY , LOAD , EXECUTE . Insert each card in numerical 
order into the card reader. After each side is loaded, “NOTE 14 ” will be displayed, with 
the exception of side three. Press EXECUTE and insert the next side. 

3. If side three is followed by a “NOTE 14 ” display, the calculator has misread 
program cards. Press ERASE MEMORY and reload the MET routine. 

4. Place graph on plotter for meteorological profile (use when reprocessing data). 

B. Data Cards 

Input data can be recorded in converted form on magnetic cards for reference to 
a particular case study. Magnetic cards may also be used to enter data into the calculator. 
These operations center around the “LOP OPT ” display and will be discussed later. 
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C. Entering Data Through the Keyboard 


1 . During processing the calculator will idle, generating a display for data entry. 

2. Following the display, type in the corresponding data and press RUN 
PROGRAM . 

Example: 

“TIME ” is displayed. 

Enter the time through the keyboard and press RUN PROGRAM . 

D. Initiating the MET Routine Loading Data by Cards 

1. Press RUN PROGRAM until “Z” is displayed. Press SIP, GTO 26 . RUN 

PROGRAM . This jumps over data input and conversion routines required for raw 
information. will be displayed. Surface density is stored on the magnetic card so 

data entry for this call is unnecessary. Press RUN PROGRAM . 

2. The card drive will activate, but in the record mode for converted raw data. 
Press RUN PROGRAM . 

3. “ LOD OPT ” will be displayed. Press LOD , EXECUTE , and load data into the 
calculator. 

4. Dry bulb temperature will not plot out with the use of data cards. 

E. Initiating the MET Routine with Raw Data 

1. Press RUN PROGRAM until is displayed. Enter the following rawinsonde 
data to the corresponding display in accordance with C, above. 

Z — Altitude (ft or m) 

0 — Wind Direction (kn or m) 

T — Temperature (dry bulb, °C) 

P — Pressure (mb) 

Each group of five parameters constitutes one data set. Fifteen data sets are 
required for processing, although nineteen sets are the maximum to be stored. In cases of 
missing data points, the following may be substituted. 
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Missing Parameter 


Substitute 


Altitude 0 

Wind Direction 1000 

Wind Speed - 1 

In the event that there are less then 15 data sets, this number can be generated by 
cycling through with RUN PROGRAM until 15 sets are generated. 

2. If data were miskeyed into the calculator, see subsection III.C for data 
correction or press STP , GTO 0 , RUN PROGRAM , and reenter data sets. 

3. After storing the rawinsonde data, press STP, GTO 12 , RUN PROGRAM . 
“METRIC? ” will be displayed. If rawinsonde data were in metric units, press GTO 19 , 
RUN PROGRAM to jump over the conversion routine; otherwise press RUN PROGRAM . 

4. Dry bulb temperature will plot out. Change pen colors. 

5. “PS” will be displayed. Enter the surface density in g/m^. 

6. The card drive will activate. Converted data can be recorded by inserting a 
data card into the card drive; otherwise, press STP , RUN PROGRAM (when initially 
using the program, it is recommended that data be recorded). A “NOTE 14 ” display 
afterwards indicates a data error or a recorder malfunction. In such case, press STP . 
EXECUTE . GTO 26 . RUN PROGRAM . With the “DS” display, press RUN PROGRAM 
and again record the data. After the second appearance of “NOTE 14 ”, if any, see 
subsection III.C for data correction. 

7. “LOD OPT ” will be displayed. This concerns only data card loading and it 
should not be associated with raw data. 

F. Plotting the Meteorology Profile 

1. After the “ LOD OPT ” display, either with keyboard or card loaded data, 
press RUN PROGRAM . 

2. Wind speed will plot out. Change pen colors and press RUN PROGRAM . 

3. Potential temperature and adiabatic potential temperature gradient will plot 
out. Change pen colors and press RUN PROGRAM . 

4. Wind direction and its scaling will plot out. 

5. “MONTH ” will be displayed. Enter in the date of the rawinsonde data with 
the following displays; “MONTH ”, “DAY”, “ YEAR” , “ TIME ”. 
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6. Date of rawinsonde data and surface temperature, pressure, and density will 
plot out. 


7. The “C-GRAPH ’’ displays the end of the meteorology profile routine. 

8. For safety, enter ^ ,D,A, and EXECUTE . This enables 

returning to this point without rerunning the routine. 


PART II. CLOUD RISE 

1. Load cloud rise program (see Part I, Section LA). There are two sides. 

2. Set up the plotter. Press RUN PROGRAM . If an adiabatic condition exists, 
plotting will commence until the 900-sec limit is reached or until a stable condition 
prevails. 

3. Once processing has stopped or if plotting did not initiate after pressing RUN 
PROGRAM at the start, change pen colors and press RUN PROGRAM until the 
completion of the plot. For each time the program stops, a numerical readout will be 
displayed, which is the time factor calculated for a given altitude. This effect indicates 
the processing of a stable condition. 

4. After completion of the printing, ‘‘Stable Atmosphere” will be printed out 
beneath the line segments which it represents. Here the program will stop for a pen color 
change. Press RUN PROGRAM to resume processing. Dates of rawinsonde data, time to 
maximum cloud rise, and “adiabatic atmosphere” (if one exists) will be plotted out. This 
completes the graph and the MET routine. 

5. If a note statement appears or the cloud rise plot seems unrealistic in 
accordance with the meteorology profile, data are in error. Most frequently trouble will 
arise in registers Rl, the surface temperature in Kelvin; R2, the surface pressure in mb; 
and R3, the surface density in g/m^. See subsection III.C for data correction and for 
recycling the program. 

6. Again, for safety, do Part I, Section LF.8. 


PART III. DIFFUSION ROUTINE 

SECTION I. CENTERLINE CONCENTRATION AND DOSAGE 

1. The cloud diffusion routine requires two cards. Three sides are used for the 
concentration and dosage routines and one side for the FA subroutines. MEMORY 
ERASE should not be depressed, since data stored in the previous routine are used to 
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begin loading. Press DBF FA , EXECUTE , LOD , EXECUTE , and insert the FA subroutine 
program card. This is stored into the FA register and may be called back for inspection 
by pressing DEF FA , LIST or DEF FA , RECALL. To load the cloud routine, press END . 
EXECUTE , and insert the cloud routine program cards. The END function separates the 
program mode from the subroutine mode. 

2. Set up the plotter for the centerline graph. Press RUN PROGRAM until “£T” 
appears in the display. Select from the meteorology profile the altitude, wind speed, and 
wind direction for the top of the surface layer. With each of the following displays, enter 
its corresponding data: 

OT — Wind direction at the top of the layer 

UT — Wind speed at the top of the layer 

OAR — Wind deviation at ground level 

HL — Height of the layer 

3. After the last data input, plotting will initiate. If, for some period of time the 
calculator should idle with no apparent output, storage data are in error. See subsection 
III.C for data correction and for recycling the program. 

4. Following the concentration and dosage plot, the height of the surface layer, 
maximum concentration, and the date of the rawinsonde data will be printed out on the 
graph. The “ ISO-GRAPH” display indicates the end of this part of the cloud routine. 

SECTION II. ISOPLETH 

1. When “ISO-GRAPH ” is displayed, set up the plotter for the isopleth graph, 
change pen colors, and press RUN PROGRAM . 

2. Cloud path will plot out in a relatively straight line, with the majority of its 
deviations centered around the origin. If there are any signs of irregularities similar to a 
spike or a jump in the line, miscalculated data were stored in memory. To continue with 
processing will lead to inaccuracies. For data correction and to recycle the program, see 
subsection III.C. 

3. In many cases the cloud path will go off scale and the arbitrary down range 
time value will be missing. This is a processing error. No action is necessary since this has 
no effect on the output. To prevent this, press STP before the cloud path reaches the 
limits of the graph, then press GTQ 25 , RUN PROGRAM . 

4. Following the cloud path plot, time of cloud rise, an arbitrary down range 
time factor, and the position of maximum concentration are marked out. The maximum 
concentration parameter and date of the rawinsonde data are then plotted out. 
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5. “C” will be displayed. Enter a concentration greater than 0.1, for 
practicability, and no less than maximum concentration. 

6. The time before initial plotting is dependent on the proximity of the entered 
value to maximum concentration. This lengthens as the two values approach each other. 
In some cases, a note statement will be displayed, indicating that the calculator is unable 
to determine the initial plotting location. Enter a new concentration value by pressing 
STP , EXECUTE , END , EXECUTE , GTO 31 , RUN PROGRAM . If only one point is 
plotted and the calculator idles for a period of time, enter a new concentration value. 

7. During plotting, if the isopleth should go off scale, the program will recycle 
for the entry of a new concentration value. If the isopleth remains within the limits of 
the graph, ‘ ‘NOTE 2 ” .will usually appear. Press RUN PROGRAM to complete the 
isopleth. If “NOTE 2 ” appears a second time, press STP, EXECUTE , GTO 31 , RUN 
PROGRAM . will be displayed. 

8. To initialize an isopleth with a concentration value, press SET/CLEAR FLAG 
during its plot. 

9. This completes the isopleth and the cloud routine. 

Subsection C affords some helpful tips on how to correct errors and shortcut data 
processing steps. 


C. Corrective Procedures and Operational Shortcuts 

Since errors can be introduced into the final results by incorrect inputs to the 
NASA/MSFC Surface Layer Diffusion Program, a careful check should be made of each 
graph for obvious errors — especially, the meteorological profile. If an error is detected, it 
is time saving to be able to correct the error and continue processing without having to 
go to the beginning and start the program at step 1. This subsection discusses the 
potential errors and how to correct them. Then the necessary steps to recycle the 
program are given. 

Since there may be times when only the output of a specific section is desired, 
the methods for bypassing unwanted output are also considered. Appendix B gives the 
location of the meteorological, cloud rise, and diffusion algorithms used in the program 
along with a flow chart and listing of the program. 

The need for a data correction can be indicated by the appearance of a note 
statement (not normal to the operation) or the results of an obviously questionable 
output which are the indications of erroneous data in storage. Mistakes in keying data 
into the calculator are generally the source of the error, although on rare occasions a 
misread in the program card would lead to a miscalculation. The case run can be salvaged 
if the error is located and corrected. 
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• With the display of “ NOTE A in B ”, where A is the note type and B is the 
line number where the error takes place, one should identify the note and its possible 
cause. 


NOTE 1: 
NOTE 2: 
NOTE 3: 
NOTE 4; 
NOTE 5: 

register. 

NOTE 6: 

register. 

NOTE 8: 
NOTE 9: 
NOTE 10: 
NOTE 12: 
NOTE 14: 
NOTE 16: 


Illegal instruction. Program card misread. 

Improper store or taking the square root of a negative number. 

Illegal instruction. Program card misread. 

Illegal instruction. Program card misread. 

Miscount in data increment calling for a nonexisting or program 

Miscount in data increment storing a value in a nonexisting or program 

Recalling nonexisting program line. 

Illegal instruction. Program card misread. 

Underflow or overflow. 

Illegal instruction. Program card misread. 

Magnetic card operation incomplete. 

Printer out of paper. 


• To correct program card misreads, press LIST , following the note statement. 
Inspect the printout with the listings given in Appendix B. Once the invalid line has been 
located, press GTO B (B is the program line), RECALL , BACK , and press DELETE until 
only the line number remains (ir. “64:”). Type in the corresponding command from 
Appendix A and press STORE (see item 1 , Recycling Meteorology Profile, to recycle the 
routine). If the program card misread is extensive, press ERASE MEMORY and restart 
the program. 

• For a single invalid data point, press GTO B , where B is the line number from 
the note statement, and LIST . Allow a few lines to print out and press STP. 

• The first line in the listing contains some register (RO, R1 . . . R103, A, B, C, 
X, Y, Z) and the operation where the error developed. The content of these registers can 
be determined, for example, by pressing B, EXECUTE , where B is the register of interest. 
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B. 


• To correct data, press N-»B , EXECUTE , where N is the new value for register 


• If invalid data cannot be located by the above method or there are a number 
of data points at fault, a data dump will be required. To initiate, press CH-X , EXECUTE , 
PRT RX : JMP (X+l^X) 100 , and press EXECUTE until approximately 100 figures have 
been printed out. From RO down to R103 will be listed consecutively, and in this way 
irregularities in groups of data may be found. To insert a correction, press N-»B , 
EXECUTE where N is the new value for register B. 

• Following data correction, the program should be recycled to the beginning of 
the routine (see Part I, Section I.E, step 3 of the program). In cases where data 
correction cannot be handled, press ERASE MEMORY and restart the program. 

Portions of the program, where an output is not needed, may be jumped during 
processing if the data required are accessible. The meteorology profile, cloud rise, and 
isopleth are a waste of operating time if the only output required is the centerline 
concentration and dosage. The method described below will show how this may be done. 

1 . Cloud Rise Plot 

a. To plot cloud rise, the MET routine should be loaded into the calculator and 
the converted rawinsonde data must be present. If data are absent from storage, they 
should be entered by magnetic card at the “LQD OPT ”, displayed or keyed in at the 
beginning of the MET routine, and followed through until “LQD OPT ” appears. 

b. After the “ LOP OPT ” display, press GTO 58 , RUN PROGRAM. Cloud rise, 
will plot out, but the date of rawinsonde data will be in error, unless Month R90, Day 
R91, Year R92, Time R93 are entered. 

2. Centerline Concentration and Dosage 

a. Load the cloud routine in the following sequence: 

Press MEMORY ERASE 

Press DEF FA , EXECUTE , LOP , EXECUTE. Insert FA subroutine. 

Press END , EXECUTE , LOP , EXECUTE . Insert cloud routine. 

b. Rawindsonde data are not required, but the following data are, although they 
are not called for by the program. 
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Data 


Register 


Surface Temperature, °K 

R1 

Surface Pressure, mb 

R2 

Time of Maximum Cloud Rise 

R85 

Maximum Cloud Height 

R86 


Enter uncalled data by pressing N-»-RB . EXECUTE where N is the data and B is 
the corresponding register number. 

c. Press RUN PROGRAM until "‘dT” is displayed. Continue operations with 
instructions in Part III, Section II of the program. 

d. Date of rawinsonde data will be in error. 

3. Isopleth 

a. Isopleth can only be generated by running the entire program. Data 
requirements are too large to calculate and enter. There are, however, methods to reduce 
processing time. 

b. Load in the MET routine and input data by card or keyboard. Follow the 
meteorology profile instructions, Part I of the program, until the “LOP OPT ” stage ends. 
Press GTO 58 , RUN PROGRAM . 

c. Let the cloud rise routine run out completely. 

d. Follow through with instructions in Part III, Section I of the program for 
centerline concentrations and dosage. 

e. After the plot has passed the point of maximum concentration, press END , 
EXECUTE , STP , EXECUTE , GTO 13 , RUN PROGRAM . 

f. Follow the instructions for isopleth, Part III, Section II of the program. 

Recycling sections of the program, following a data correction, is required to reset 
calculations with the altered data values. Simple jumps from the point of correction to 
the start of a routine eliminate the need to run several sections of the program to 
generate an output. The commands listed below may be used in any point of a given 
routine for recycling back to its beginning unless stated otherwise. 
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1 . Recycling Meteorology Profile 

a. While entering raw data, if a miskeyed data point was input, two methods 
may be taken for correction. Either search for the data point and correct it, as in Part I, 
Section I.E, of the program and continue to enter data sets, or wipe out stored data and 
reenter the data sets by pressing STP, GTO 0 , RUN PROGRAM. 

b. If more than 19 data sets were entered, command lines are erased, rendering 
the program and data unsalvageable. Press ERASE MEMORY and reload the MET 
routine. 


c. To recycle during the plotting sequence, press STP, 4-^X , EXECUTE , GTO 27 , 
RUN PROGRAM . If the graph is to be replotted, dry bulb temperature will be missing. 

2. Recychng Cloud Rise 

Processing should be in the stop mode. Press R1-273.15*^R1, EXECUTE, GTO 
RUN PROGRAM . ~ 

3. Recycling Centerline Concentration and Dosage 

Press EXECUTE , END , EXECUTE, GTO 0 , if new top of the layer 

parameters are to be entered, otherwise GTO 2 , RUN PROGRAM . 

4. Recycling Isopleth 

a. In the cloud path mode, press R97^R64, EXECUTE, GTO 13, RUN 
PROGRAM . 

b. In the isopleth mode, press STP, EXECUTE , END , EXECUTE , GTO 13 , RUN 
PROGRAM . 

5. General recycling can be accomplished if at the completion of a subroutine 
the data are recorded. This is accomplished by REC ,”,D,A,'\\ R105 , and EXECUTE . 
This permits the investigator to return to this point and reload the data by pressing 
LOAD , EXECUTE . Then the program can be rerun. This procedure is helpful in initial 
operation of this program or when doing extended investigations with this program. 

Although these corrective procedures cover most cases that will be encountered, 
there will be times when the operator of the program will want to employ his own 
techniques. While tliese procedures may normally be satisfactory, great care should be 
exercised because of the complexity of the computational interfaces. 
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SECTION IV. CONCLUSION 


The NASA/MSFC Surface Layer Diffusion Calculator Program affords a real-time, 
online transport description for the environmental dispersion of the Delta-Thor rocket 
exhaust effluents utilizing elhpsoidal dispersion in a Eulerian reference frame, as given by 
model 3 of the NASA/MSFC Multilayer Diffusion Model. This calculator program has 
been successfully employed at Cape Kennedy to make real-time predictions of the 
transport of Delta-Thor effluents for the deployment and evaluation of the Langley 
Research Center and Kennedy Space Center effluent monitoring network. 

The primary difference between the NASA/MSFC Surface Layer Diffusion 
Calculator Program and the NASA/MSFC Multilayer Diffusion Computer Program is that 
the calculator program has been specialized for a specific vehicle, the Delta-Thor, and a 
specific altitude, the surface layer. This program is similar to the one for the Titan III 
[ 1 ] . The calculator program incorporates the meteorological and cloud rise algorithms 
into one program, permitting a rawinsonde sounding to be the input; whereas, the 
computer program utilizes the results of these algorithms. The predictions obtained from 
both programs are identical for the surface effects from the effluents. The data 
preparation and reduction time for the calculator program is approximately one-half 
hour, as compared to 2 hours for the computer program. The computer program, 
however, does afford a more detailed prediction than can be obtained with the calculator 
program. 

The NASA/MSFC Surface Layer Calculator Program has features that make it a 
desirable vehicle for use in making diffusion predictions for the Delta-Thor launches, 
especially where the graphical output satisfies the prediction required. The NASA/MSFC 
Multilayer Diffusion Computer Program is a better vehicle for use in general analysis work, 
where flexibility and detail in the results are required. Thus, while the application of the 
calculator program to diffusion analysis is a step forward in online, real-time predictions, 
it does not eliminate the need for the computer program in diffusion analysis. 


38 



APPENDIX A 


SUMMARY OF THE INSTRUCTIONS FOR THE OPERATION OF THE 
NASA/MSFC SURFACE LAYER DIFFUSION CALCULATOR PROGRAM 


The objective of this summary is to provide a guide for tlie operation of the 
calculator program. The details for the initial operation of this program are given in 
Section III.B. Since, in normal operations of this calculator program, only a quick 
reminder of some key operations is all that the operator will require, the following 
summary is provided. Even this may be too detailed for many and, therefore, a flow 
chart for the operation of this calculator program has also been provided (Fig. A-1). 

The user is reminded that the procedures for treating program errors and 
eliminating unwanted parts of the program are given in Section III.C. This discussion does 
not include the points covered in that section. 


Operation of Program 


Data called for by the machine are entered in the following sequence. 


“DISPLAY” The program makes the call for data, type in data, and press RUN 
PROGRAM . 

MET Profile PART I. MET ROUTINE 


Load Program 


1. This routine consists of two cards, three sides. Press 
MEMORY ERASE and load in the routine. 


Loading Data 2a. If data are to be entered by magnetic cards, press RUN 

by Card P ROGRAM until “Z” appears in the display. Press STP , GTO 

26 , RUN PROGRAM . “PS” will appear in tlie display though 
data entry for it is not necessary. Press RUN PROGRAM . The 
card reader will activate in the record mode for data records, 
which at this point do not exist. Press STP , RUN PROGRAM . 
“LOD OPT ” will be displayed. Press LOD , EXECUTE , and 
insert data card. Dry bulb temperature will not be plotted 
with data card usage. Press RUN PROGRAM to resume 
processing. 


Entering Raw Data 2b. The met routine should be loaded. Press RUN PROGRAM 

until “Z” appears in the display. Enter data with the 
corresponding display. “Z” altitude, “0” wind direction, “V” 
wind speed, “I” temperature, and “P” pressure. No more than 
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Figure A-1. Flow diagram for program operations. 













Recording Data 


Cloud Rise 
Load Program 


95 values may be entered. For missing data, the following 
may be substituted: Altitude, 0; wind direction, 1000; wind 
speed, -1. 


3. After entering the last data value, press STP , GTO 12 ,. RUN 
PROGRAM . “METRIC? ” whl be displayed. If data were in 
metric units, press GTO 19, RUN PROGRAM; otherwise press 
RUN PROGRAM . 

4. Dry bulb temperature will plot out. “DS ” will appear in the 
display. Set up plotter. Enter surface density value, change 
pen color, and press “RUN PROGRAM ”. 

5. The card drive will activate. If converted data are to be 
recorded, insert a magnetic card; otherwise press STP , RUN 
PROGRAM . 

6. “LOP OPT ” will be displayed. Press RUN PROGRAM. 

7. Wind speed will be plotted out, stopping afterwards for a pen 
color change. Press RUN PROGRAM . 

8. Potential temperature will be plotted out, stopping afterwards 
for a pen color change. Press RUN PROGRAM . 

9. Wind direction, wind direction scale, surface conditions, and 
date of meteorology data will be plotted out. “C-GRAPH ” 
will be displayed. 

PART II. CLOUD RISE ROUTINE 

1. This routine consists of one card, two sides. Load routine. 

2. Center graph paper on plotter and press RUN PROGRAM . 

3. During processing, if the plotter should stop, change pen color 
and press RUN PROGRAM as many times as needed to 
complete the graph. 

4. After the plotting of “STABLE ATMOSPHERE ”, change pen 
color and press RUN PROGRAM . This completes the MET 
routine. 
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Centerline Concen- 
tration and Dosage 
Loading Cloud 
Routine 


Isopleth 


PART III. CENTERLINE CONCENTRATION AND DOSAGE 
CLOUD ROUTINE 

1. Load cloud routine, two cards, three sides, and FA 
subroutine, one card, one side, in the following sequence. Do 
not press MEMORY ERASE ; press END, EXECUTE . LOP , 
EXECUTE . Load cloud routine. 

2. Press RUN PROGRAM until “OT” is displayed. Enter the 
corresponding data to each display: OT — wind direction at 
the top of the layer, UT — wind speed at the top of the 
layer, HL — height of layer, OAR — wind deviation at 
reference. 

3. After the final data entry, the program will plot concentration 
and dosage and then print- out height of layer, maximum 
concentration, and date of meteorology data. This completes 
the concentration and dosage graph. An “ ISO-GRAPH ” will 
be displayed. 

PART IV. ISOPLETH CLOUD ROUTINE 

1. Center graph paper on plotter and press cloud routine RUN 
PROGRAM . 

2. Qoud path, maximum concentration parameters, and date of 
meteorology data will be plotted out. 

3. “C^’ will appear in the display. Enter in a concentration value 
greater than 0.1 but less than maximum concentration. The 
corresponding isopleth will be plotted. 

4. To initialize the isopleth during plotting with a concentration 
value, press SET/CLEAR FLAG . 

5. If “NOTE 2” is displayed, press RUN PROGRAM to 
complete the plot. After the second appearance of “ NOTE 
1”, press GTO 31 , RUN PROGRAM . “C” will appear in the 
display. 

6. If “NOTE 2 ” did not appear, the program cycles back to plot 
another isopleth, displaying “Q” for a new concentration 
value. 

7. This completes the cloud routine. For completion of each 
graph there is a corresponding routine for axis and lettering. 
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APPENDIX B 


THE NASA/MSFC SURFACE LAYER DIFFUSION CALCULATOR PROGRAM 

LISTING AND ALGORITHMS 


A listing of the NASA/MSFC Surface Layer Diffusion Calculator Program, along 
with flow charts, is given in this appendix. In addition, the location of the algorithms 
used in the program is presented. 

This calculator program was designed for use on the HP9820 desk calculator, with 
option 001 (429 total registers). The user definable functions I-III ROM are required in 
this program. The program requires all the storage available in this configuration. A listing 
of the MET routine is given in Table B-1, the cloud rise routine is given in Table B-2, the 
diffusion routine is given in Table B-3, and the FA subroutine is given in Table B-4. 
Figures B-1 through B-8 are the flow diagrams for the program. 

Since experience in making diffusion predictions generally leads to a desire to 
update the algorithms utilized, the meteorological, cloud rise, and diffusion algorithms are 
given in Table B-5 and B-6. However, the algorithms employed in this program have 
evolved from years of empirical experience in monitoring the dispersion of the Delta-Thor 
effluents and, therefore, should require practically no changes. 

To assist the user in making his initial predictions with this program and to 
permit the checking of its operation, some average values have been included for the 
parameters in Table B-7. Table B-8 provides the four programs for formating and labeUng 
each of the graphs. In general, time can be saved by making a master for each graph and 
reproducing it for use with the diffusion program. 
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TABLE B-L LISTING OF MET ROUTINE 
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TABLE B-2. LISTING OF CLOUD RISE ROUTINE 
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TABLE B-3. DIFFUSION ROUTINE 
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TABLE B4. FA SUBROUTINE 
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ENTER DATA SET 


FLG ^ ^ 

SET FLG 

|YES 





SHIFT REGISTORS 
FROM SET (A) UP 


TEST SET (A+l^A 


INSERT SET (N) 


L = LINE NUMBER 
SET (N) = THE ENTERED DATA SET 
SET (N-1) = LAST DATA SET STORED 

SET (A) = DATA SET, INITIATING AT THE FIRST DATA SET 


Figure B-1. Flow diagram for sort subroutine. 






LI 3 



Figure B-2. Flow diagram for convert subroutine. 
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127,28,29,30 


1 


PLOT 

WIND SPEED AND LABELING 


L31,32 


I 


PLOT POTENTIAL TEMP. 
AND LABELING 


L33,34,35 


I 


PLOT ADIABATIC POTENTIAL 
TEMPERATURE GRADIENT 


L36,37,38,39,40 


I 


CALCULATE MINIMUM, 
MAXIMUM WIND DIRECTION 


L41 ,42,43,44,45 


I 


CALCULATE WIND 
DIRECTION SCALING 


L46 


I 


PLOT WIND DIRECTION 
SCALING 


L47,48,49,50,51 


I 


PLOT WIND DIRECTION 


L52 


I 


PLOT SURFACE 
TEMPERATURE 


L53 


I 


PLOT SURFACE 
PRESSURE 


L54 


I 


ENTER RAWINSONDE 
DATE 


L55,56 


I 


PLOT RAWINSONDE DATE 
(ON ALL GRAPHS) 


L57 


I 


PLOT SURFACE 
DENSITY 


I 


END I 


L = LINE NUMBER 


Figure B-3. Flow diagram for meteorology profile subroutine. 



END II 


Figure B-4. Flow diagram for cloud rise routine. 
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Figure B-5. Flow diagram for centerline concentration and dosage routine 












4 


'•'V--1«Sing\»es 


INTERPOLATE 
WIND p E EP 

U6 I 

INTERPOLATE 
WIND DIRECTION 


CALCULATE 
DOWN RANGE 
DISTANCE 


CALCULATE 

DIRECTION 


L^l T 

CONVERT 
POLAR TO 
RECTANGULA R 

L21 i 


CLOUD NO 

stabilization'>^ 

J^EACHED 


NEXT DATA 
SET 


L23 T 
STORE CLOUD 
STABILIZATION 
COORDINATE _ 

L24___|[ 

CALC. DISTANCE 
AFTER 

STABILIZATION 


CALC. DIRECTION 
AFTER 

STABILIZATION 


PLOT MARKER. 
TIME 


L27 

/ PLOT MARKER 
f FOR MAXIMUM 
CONCENTRATION 


L = LINE NUMBER 


L2B T 
/ PLOT TIME OF 
I MAXIMUM 

CONCENTRATION 


L29 ? 

/plot DIsf^CE AND- 
/angle TO maximum/ 

f CONCENTRATION / 


Figure B-6. Flow diagram for cloud path routine. 
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ITER 



LINE NUMBER 


Figure B-7. Flow diagram for isopleth routine. 









LINE NUMBER 


Figure B-8. Flow diagram for FA subroutine, 

















TABLE B-7. AVERAGE VALUES FOR SURFACE LAYER DIFFUSION ROUTINE 


‘^yo’ ‘^xo’ °zo 

194.7 




m 

-0.0979 




^A 

-0.07203 




Down range distance 

0 

500 

1000 

1500 

a^Jz 

7.58 • 10^ 

1.056 • 10® 

1.430 • 10® 

V) 

o 

00 

S 

0 

2 • 10'® 

4 • 10'^ 

2.61 • lO"- 

Vertical term 

0.00163 

0.01215 

0.04329 

0.10263 



TABLE B-8. LABELING ROUTINES FOR GRAPHS 


0 : 
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^^400-h RXE 0? 0? 1 ? 
1 0 0 5 H /: E 0 ? 2 D 0 3 1 
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